A t her mo dynamical approach to dissipation 

range turbulence 

Jacopo Bellazzini 



Dipartimento di Ingegneria Aerospaziale 
Universita di Pisa 
Via G. Caruso 56100 PISA -Italy 

Abstract 

A model to explain the statistics of the velocity gradients in the 
dissipation range of a turbulent flow is presented. The experimentally 
observed non-gaussian statistics is theoretically predicted by means of 
a thermodynamical analogy using the maximum entropy principle of 
ordinary statistical mechanics. 

PACS numbers: 47.27. Gs 

The dynamics of fluids at high Reynolds number is one of the most interest- 
ing subject of statistical physics. In the limit of inflnite Reynolds number, 
the celebrated theory of Kolmogorov (K41) IJ suggests that the small scales 
statistics is characterized only by the mean rate of energy dissipation per unit 
mass e and the scale /. The K41 theory is based upon the concept of self- 
similarity of the inertial range [Id, L] and implies that the velocity gradients 
Svi{x) = \v{x + I) — v{x)\ scale as 

< 6vf >~ (e/)^/l (1) 

The experimental evidence of the breakdown of eq.(IT)), see 013111, due 
to the presence of intermittency, induced Kolmogorov and Obukhov IH] to 
modify the K41 theory introducing the fluctuation of the mean flux of energy 
e. Indeed, as pointed out by Siggia [7J, the turbulent flow can be described as 
the superposition of two fluid states, the coherent structures and the random 
fluctuations. Hence, in the inertial range, the flow preserves organized struc- 
tures as, for instance, vortex sheets and fllaments. The coherent structures 
are "rare" events that contribute to the non-gaussian high-amplitude fluctua- 
tions. Many models for such superposition have been proposed, either based 
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on a mapping-closure theory or on the fractal approach jHJ 13 UHl IHl 1121 CHI ■ 

However, also at very small scales I ~ l^, when the coherent structures are 
annihilated by dissipation, the tails of the velocity gradients statistics are 
deeply non gaussian. Indeed, experimental results show that the tails of 
the probability density function (PDF) of the velocity gradients are close to 
exponential |T^ for small /. In this paper, we conjecture that the PDF of 
the velocity gradients in the dissipation range, / ~ Z^, can be described by 
means of ordinary statistical mechanics. We propose to consider the velocity 
gradients Svi of the random incoherent flow like the momentum p of a par- 
ticle in a ideal gas with a "temperature" T that changes slowly. Specifically, 
we consider this temperature not to be fixed but x^-distributed with degree 
n = 3; the resulting PDF agrees both with the experimental results and with 
the theoretical structure function. 

In the ordinary description of an ideal gas, the momentum p of a particle 
with mass m may be considered as a random variable. From the maximum 
entropy principle one obtains the density distribution 



where Z is the partition function and P = Here T is the temperature of 
the thermal bath and k is Boltzmann's constant. 

Now, we consider the temperature of the termal bath not to be constant in 
time, but fluctuating with a given distribution f(T). In practice, we suppose 
that the time scale on which T fluctuates is much larger then the typical 
time to reach equilibrium. Let V^(p|T) be the conditional probability, and let 
■^(p) be the probability to observe a certain value of p no matter what T is. 
Then, the following relation holds 



Due to the maximum entropy principle of ordinary statistical mechanics, the 
conditional probability ipiplT) is gaussian but the probabihty to observe a 
certain value of p is not. 

The main idea of the paper is to consider the velocity gradients Svi{x) in 
the dissipation range and in the limit of infinite Reynolds number, as the 
momentum p of a particle fulfilling eq.Q with m = k = 1. In our approach, 
the "temperature" corresponds to the averaged kinetic energy fluctuation 
< 6vf > and will be referred with the symbol T^. 

As shown by Kraichnan if we observe the statistics of the velocity gra- 
dients in the inertial range, considering both the coherent structures and the 



^(p) = -e 



(2) 




(3) 
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random fluctuations, the equipartition of energy does not apply. However, in 
the dissipation range, if we suppose that the fluctuations of the velocity gra- 
dients are given only by the incoherent fluctuations, we argue (and assume) 
that the equipartition of energy holds. The physical implication of this model 
is that the velocity fluctuations, in the dissipation range, may be described as 
a cascade of thermal baths. At every scale the ordinary statistical mechanics 
holds, i.e. the conditional probability for a fixed temperature Tt is gaussian, 
but, due to the changes of such temperature, the final distribution will be 
different from the gaussian one. 

Due to the three dimensional nature of the turbulent flow, we suggest that 
the temperature Tt fluctuates with a x^- distribution with degree n = 3. Ac- 
tually, such physical model for the distribution of the averaged kinetic energy 
fluctuations /(Tj) has been proposed recently by Beck J6j in the formalism 
of nonextensive statistical mechanics. 

The temperature Tt has to be considered as a random variable given by the 
expression 
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Tt = Y.^"^ (4) 

i=l 

where Xj are independent gaussian variables with zero average. The mean 
value of the averaged kinetic energy fluctuations Tt depends on the scale /. In 
particular, at scale /, the mean value < Tt > is 3g{l), where g{l) =< >. 
This means that the probability density /(Tj) is given by 

If we introduce eq.® into eq.® we obtain the following integral 

(J r+oo svf Tt 
iIj{6vi) = —TTTTj:: I e'^^e'^dTt, (6) 



where C = (4yvrr(3/2))~^ is the normalization constant. The latter integral 
can be evaluated exactly. Indeed, eq.® can be rewritten as 

^(M) = ^^Sv^K, (-ip^) , (7) 

where Ki{-) is the modifled Bessel function of the second kind. 
Eq-© provides a theoretical prediction of the statistics of the velocity gra- 
dients of the incoherent turbulence in the limit of inflnite Reynolds number. 
Obviously, the statistics of eq.((H) depends on the scaling law g{l) of the mean 
temperature < Tt >. 
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Given the relation between the mean temperature < Tf > and the scale /, 
the moments of the velocity gradients can be evaluated easily by means of 
eq.((Z)). Indeed, we obtain 

< 5vf >= f{q)g{iy'\ (8) 
where /(g) is a function of q depending on the moments of the Bessel function 

f{q)=2C I y^+'K^{y)dy. (9) 

Equation (jH)) shows that the moments of the velocity gradients depend only 
on the scaling law of the mean temperature < Tt >. 

At the very small scales, / ~ Id, the velocity field is smooth due to the lack of 
singularities. Therefore, the structure function scales as < 6vf >r^ l'^. Then, 
due to eq.®, we propose the following scaling law 

g{l)=l\ (10) 

With this assumption, the theoretical prediction of the velocity gradients 
statistics is given by the following equation 

^(5^0 = (^^) . (11) 

We will show that such theoretical prediction is in complete agreement with 
the experimental results. 

Actually, experimental data in the dissipation range are difficult to acquire. 
We show the comparison between the theoretical prediction and the exper- 
imental data of Van de Water These experimental data correspond to 
the transverse PDF of a turbulent flow behind a grid in a closed wind tunnel 
at a nearly-dissipative scale {j- ~ 10). Figure^shows the agreement between 
the experimental results and the theoretical prediction. 

Although the moments of the velocity gradients distribution of the inco- 
herent turbulence depend on the scaling law of the mean temperature < Tt >, 
the probabilty ip{6vi) has an asymptotic decay close to the exponential func- 
tion whatever is the relation between the scale / and the mean temperature 
< Tt >. Indeed, the asymptotic shape of the velocity gradients is given by 

ip{6vi)r^e^, (12) 

where c(/) is a function of the scale / that depends only on g{l). 

As shown by Benzi et al. the PDF of the velocity gradients observed 
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Figure 1: Log-linear plot of the probability distributions ip{Svi) of the velocity 
gradients. The theoretical prediction is obtained with 1=0.58 (dashed- dotted 
line). The experimental data of Ref. fT^ are indicated by symbols. 



experimentally in the inertial range can be obtained by a non-trivial super- 
position of stretched exponentials, corresponding to the various exponents of 
the multifractal approach. Moreover, such distribution is explicitly depen- 
dent on the Reynolds number. On the contrary, at very small scales / ~ Z^, 
the tails of the PDF are exponential-like and independent of the Reynolds 
number. In practice, at very small scales, the exponential like decay of the 
PDF is a universal feature that, using the present approach, may be described 
by means of ordinary statistical mechanics. 

In conclusion, we have presented a description of turbulence at very small 
scales / ~ /d, where the coherent structures may be assumed to have been 
annihilated by dissipation. A thermodynamical model for incoherent turbu- 
lence, providing an exact expression for the PDF (in terms of the modified 
Bessel function of the second kind), has been presented. We have shown that 
the velocity gradients fluctuations at very small scales may be described as 
the momentum p of a particle in a gas with a temperature T that fluctuates 
slowly. Hence, the present model explains how even the random incoherent 
fluctuations in the dissipative range may give rise to non-gaussian statistics. 
We hope that this effort will results in a deeper understanding of the physics 
underlying the phenomenon of fully developed turbulence and, more in gen- 
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eral, of complex processes characterized by the coexistence of equihbrium 
and non-equihbrium conditions. 



The author is very grcatful to G. Buresti for fruitful discussions and the 
continuous help in the preparation of the paper. 
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